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Abstract. In a recent paper Coldea et al (2010 Science 327 177) report observation 
of the weak confinement of kinks in the Ising spin chain ferromagnet CoNb206 at 
low temperatures. To interpret the entire spectra of magnetic excitations measured 
via neutron scattering, they introduce a phenomenological model, which takes into 
account only the two-kink configurations of the spin chain. We present the exact 
solution of this model. The explicit expressions for the two-kink bound-state energy 
spectra and for the relative intensities of neutron scattering on these magnetic modes 
are obtained in terms of the Bessel function. 
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1. Introduction 

Very recently Coldea et al [1] reported the impressive results of inelastic neutron 
scattering experiments on the quasi-lD ferromagnetic CoNb206 (cobalt niobate) single 
crystal. The essential physics of this material at low temperatures can be described in 
terms of the quantum Ising spin chain model, which is paradigmatic for the theory of 
the quantum phase transitions [2]. In presence of the magnetic field h± transverse to the 
easy magnetization axis, the ground state state of the model can be either ferromagnetic 
or paramagnetic depending on the strength of the external magnetic field h±. The 
transition between the two phases occurs (at zero temperature) at the critical value of 
the transverse field h± = he- This critical point belongs to the 2D Ising universality 
class. 

The spectra observed by Coldea et al [1] display certain very subtle features 
providing experimental confirmation of two long-standing theoretical predictions [3l H] , 
which relate to the Ising model. Directly at the critical transverse field h± = = 5.5T, 
the ratio of two lightest quasiparticles approaches to the 'golden ratio' indicating the 
hidden Eg symmetry in the critical Ising model in the longitudinal magnetic field, as it 
was predicted by A.B. Zamolodchikov |1] nearly two decades ago. On the other hand, 
at zero magnetic field, the observed energies of five lowest magnetic excitations were 
proportional to the absolute values Zn of zeroes of the Airy function, Ai(— z„) = 0, in 
agreement with the theory of the kink confinement originating in 1978 from the work 
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of McCoy and Wu [3]. On the resent developments in this field see El El El El HDl E] , 
further references can be found in the monograph [T2] . 

Confinement of topological excitations typically takes place in two dimensions (one 
spacial and one time dimension), if the discrete vacuum degeneracy is explicitly broken 
by a small interaction term. In the simplest heuristic approach, two confined kinks 
in the ferromagnetic Ising chain are treated as two quantum particles moving in a 
line — oo < x < oo and attracting one another with a linear potential A|a;|. The 
latter can be induced in the quasi- ID Ising ferromagnet either by a weak external 
longitudinal magnetic field, or by the weak coupling between the magnetic chains in 
the 3D magnetically ordered phase [H [131 IH] ■ In this approach, the relative motion of 
two kinks is described by the Schrodinger equation 

1 <f 

-f^lljn{.x) + X\x\lljn{x) = 5En1p{x), (l) 

with a skew-symmetrij^l wave function il){x) = —ip{—x). It immediately leads to the 
energy levels of the kink bound-states [3] 

6En = zj?l^m-^l'\ n = 1, 2 . . . (2) 

The simple theory of confinement based on ([2]) implies the quadratic dispersion 
law p^/ (2m) for a free kink, and ignores discreteness of the spin chain. Though these 
approximations are reasonable for small enough momenta of the composite two-kink 
bound states, a more systematic approach is required to describe their spectra in the 
whole Brillouin zone. 

The effect of the lattice discreteness on the kink confinement in the non-critical 
Ising spin chain has been studied in ref. [16] in the Bethe-Salpeter equation approach 
[71 [8]. To interpret the full experimental spectra in CoNb206, Coldea and his colleagues 
proposed a different phenomenological model, which takes into account only the two- 
kink (i.e. one-domain) configurations of the spin chain. The Hamiltonian of this model 
is given in [T7] by equation (SI), which we reproduce in equation ([3]) below. The 
neutron scattering spectra were compared by authors of ref. [1] with an approximate 
perturbative solution [18] of their phenomenological model. It is interesting, that this 
model admits an exact solution, which we describe in the present paper. Our main 
results are equations (j24l) . ( l25l) . and ( l43l) . which express the energy spectra and relative 
neutron-scattering intensities of the magnetic excitation modes for model ([3]) in terms 
of the Bessel function. 

The rest of the paper is organized as follows. In Section |2] the phenomenological 
model introduced in [17] is described and its exact energy spectrum is obtained. These 
exact spectra are analyzed in several asymptotical regimes in the small magnetic field 
limit in Section [31 Section [H contains calculation of the dynamical correlation function, 
which is proportional to the neutron-scattering intensities. The details of calculations 
are described in two Appendices. Concluding remarks are presented in Section [5l 

% Equation ([1]) with a symmetric wave function '^(x) = tp{—x) can describe the two-kink bound states 
in the 3-state Potts field theory [TSl [TO] ■ 



On the kink confinement in CoNb206 



3 



2. The two-kink model and its energy spectrum 

In this Section we study the eigenvalue problem for the model Hamiltonian defined by 
equation (SI) of [T7] : 

H\j, I) = J\j, I) - a[\j, / + 1) + / - 1) + \j + 1,1- 1))(1 - 5^,0 

+ ij - 1, / + 1)] + Mb, - -1,1) + \j + 1,1)) + p% i)^/,i- (3) 

Here \j,l) denotes the two-kink state of the ferromagnetic spin-1/2 chain, 

ij, /) = i...nm...mt...), 

the indices j and / give the starting position and the length of the down-spin cluster, 
j = 0,±1,±2,..., and / = 1,2,... Parameter J characterizes the energy needed to 
create two kinks. The terms proportional to a describe the nearest neighbour hoppings 
of kinks along the chain. The long-range attraction between the kinks is represented in 
([3]) by the term hj, where is the effective longitudinal magnetic field. The short- 
range interaction /3- and /3'-terms were introduced in [T7j to describe the experimentally 
observed 'kinetic mode' - the well localized bound-state mode near the Brillouin zone 
boundary [l]. Note, that /3- and /3'-terms has no analogue in the standard Ising spin 
chain Hamiltonian. 

In the momentum basis 

oo 

\P,l)= J2 exp(iPj)|j,0, (4) 

j=-oo 

the Hamiltonian is diagonal in the momentum variable P and acts on the basis state as 
follows 

H\P, I) = -a[(l + e^^)|P, / + 1) + (1 - 5i,,)il + e-'^)|P, /-!)] + 
[J + hJ + (/?' - 2/3 cos P)6i,i]\P, I) . (5) 
The eigenvalue problem 

HMP)) = E{P)mP)) (6) 
takes in basis (jlj) the explicit form 

[J + hJ + (/?' - 2/3 cos P)6i,i - E{P)]ij{l, P) - 

-2«cos(P/2)[V^(/ + l,P) + (l-5i,i)^(/-l,P)] = 0, (7) 

where / = 1, 2, . . ., the momentum P has the Brillouin zone — vr < P < tt, and 

27r 5{P' - P) ipil, P) = exp(-iP//2)(P', /|<I>(P)). 

Problem ([6]) can be easily solved in the case of zero magnetic field. The energy spectrum 
aX hz = has the continuous part 

E{k, P) = e{k, P) = J ~Aa cos(P/2) cos fc, (8) 

and one localized bound state, the 'kinetic mode', 

E,UP)-JHP'-2M + f^^. (9) 
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Note, that the wave functions corresponding to the continuous spectrum and to the 
kinetic bound state mode read, respectively, as 

ip{l,P;k) = Asm[k{l -I) +-f], (10) 

lPuin{lP)=Bz\ (11) 

where A and B are the normahzing constants, and 
2«cos(P/2) 



2/3cosP-/3' ' 

/3'-2/3cosP 



cot 7 = cot k + 



(12) 



2asm.k cos(P/2) 

Requirement |z| < 1 imphes, that the kinetic mode (|TT1) exists only in the region 
near the Brillouin zone boundary, at 



cos(P/2)| < 



-a + ia^ + 



+ 4/3/3')'/' 



2/3 



(13) 



E(P) 
J+4a 



J-4a 




n 
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Figure 1. The energy spectrum of ([7]) at /i^ = 0, according to 
corresponds to tlie continuous spectrum. 



2n 



The filled area 



At the edges of this region, the gap between the continuous spectrum and the kinetic 
mode vanishes. The resulting energy spectrum at /i^ = is shown in Figure [H 

Returning to the original eigenvaule problem (JTj) with /i^ > 0, we rewrite it as 
follows 



0, 



(14) 
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where / = 1, 2, . . ., the eigenfunction P) vanishes at / — > +oo, and 

E{P)~J ^ 13' -2{3cosP K . . 

4acos(P/2)' 2acos(P/2) ' ^ 4acos(P/2)" ^ ^ 

It is possible to extend equation ( IT4l) to all integer / G Z, not necessary positive. 
Let us continue ^/'(P, /) skew-symmetrically to negative / denoting 

r^(/,P), forZ = l,2,... , 

^(/,P) = I 0, for Z = 0, (16) 

i-V^(-/,P,), for/ = -1,-2,... 

One can easily check then, that the odd (in /) function \E'(Z,P) solves equation 

for I = 0, ±1, ±2, . . ., iff the function ip{P,l) is the solution of the equation (fT4l) for 
/ = 1,2,... 

After the Fourier transform, equation f lT7|) takes the form 



a / 1\ /■ dz' (/)f2'^ 

'51 



with e{z) = ~l(z + -] , (19) 



vl/(l) = (h P-U^) [--z]. (20) 



2 V z 
z ( 

. —<\>n{z) 

Here 5*1 denotes the unit circle in the complex plane, the both z and z' variables lie in 
this circle, \z\ = \z'\ = 1. Integration in f|T8|) is taken in the counter-clockwise direction 
and understood in the sense of the Cauchy principal value. The function (j){z) is defined 
as 

oo oo 

<Piz) ^Y.'' ^) = - ^"') (21) 

/=— oo 1=1 

and satisfies the symmetry property (j){l/z) = —(f){z). We have dropped the explicit 
indication of the P-dependence in (f){z) and ^(1), the full notation for these quantities 
should be (f){z,P) and ^(1,P). 

At a = 0, problem (fT8|) reduces to the First Toy Model, solved in Subsection 6.1 
of [16]. At a 7^ 0, problem (fTSjl can be solved by the same method. The result for the 
eigenvalues A„ reads as 

A„ = -/ii/„, (22) 

where Un are the solutions of the equation 

J,„(l//i) + aJ,„+i(l//i) = 0, (23) 



with J,y{x) being the Bessel function of order u. In Appendix A we describe an 
alternative proof of this result. 
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Accordingly, the solution of the eigenvalue problem (E]) reads as 

E^{P) = J-hiyn, (24) 
and z/„ are the solutions of equation 

2acos(P/2)J^„[4/i;iacos(P/2)] + (/?' - 2/3 cosP) J^„+i[4/i;^a cos(P/2)] = 0. (25) 

Figure |2] shows 30 lowest modes, calculated from fj2^ . ( 125|) with the Hamiltonian 
parameter values 

J = 1.94 meV, a = 0.12 J, /3 = 0.17 J, /3' = 0.21 J, = 0.02 J, (26) 
chosen by Coldea et al [17] to give the best fit of the experimental results. 



E„iP) (meV) 



3.0 



2.5 



2.0 



1.5 



••••••••••••••••••••••••••••••••••••••• 



(•••••••••• 
•••••••••• 

•••••••••••• 

i ••:!! •••j: 



•••• 

f 

^ 



«•••••••••• 

•••••••••* 



i!!! 



••• 




•• ••••••• 



,«••• 













iiiii:::: •••••••••••••••••••••••••••••••••••••••••••• 









• •••••••••••{••••••••••••••••••••••Ilaaaatmtt 



••••••••• 




••••• 















••. •.. 

XS:::>:::: 



••••••• 

••••••• 



••••••• 



TT 

2 



3_n 
4 



T ^ - 27r 



In 
2 



4 



Figure 2. Energy spectra of 30 lightest niodes calculated from (|24|) . (|25|) . (|26 



3. Weak-Zi^ expansions 

If the longitudinal magnetic field is weak hz — ?■ +0, the long range confining potential 
hzl between the kinks becomes small. Corresponding asymptotic expansions for the 
spectra En{P) of their bound sates can be extracted from the obtained exact results, 
either by means of appropriate asymptotic formulas for the Bessel functions, or by direct 
calculation of the integral in the left-hand side of equation [see equation (]A.11|) 

l^dz' (^l + a^exp|^[^(y)]|=0, (27) 
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by the saddle point method at /i — )■ +0. Here J-'{z) = — Alog^ — ^{z — z and 
contour C is shown in Fig. lAll in [Appendix "A| In this hmit, one can distinguish several 



asymptotical regimes for different n and P, depending on the location of the saddle 
points 

zl, = -A„, ± ^/Xl-l (28) 

in the 2;'-plane in integral (!27|) . Since these calculations are very similar to those 
described in [16], we present here only the final results. 

3.1. Low energy expansion 

If ~ 1 and P is well below the Brillouin zone boundary, A becomes close to —1, and the 
saddle points fl28|) merge at z[ 2 = 1- ^^is case, one can use the low energy expansion 
in fractional powers of hz'. 

E^{P) = J-Aa cos(P/2) + [2a cos{P/2)Y'^hl'^Zn + (29) 

K{P'-2PcosP) hTzl I 

/3'-2/3cosP + 2acos(P/2) 60 [2a cos(P/2)]i/3 ^' 

with {—Zn) being the zeros of the Airy function, Ai(— 2;„) = 0. 
The analogous expansion for A„ reads as 

+ ^ + + (30) 

The two leading term in expansion (129]) can be written in the form 

92 e{k,py^'^ 



K(P) = £(0,P) + 



hl'^Zn + OiK), (31) 



fc=0 

where e{k,P) is the continuous two-particle spectrum at = given by ([8]). This 
formula being in agreement with equation (87) of ref. [12], is to a large extent 'model 
independent'. That is, its applicability does not depend on the explicit form of the 
two-particle spectrum e{k,P), provided that the factor in the square brackets in ( l3Ti) 
is positive. Note, that the second term in ( 13T|) explicitly depends on the bound-state 
momentum P, in contrast to the finite momentum formula proposed by Coldea et al to 
generalize equation (3) of ref. [1], see the in line equation for mj{k) in Page 8 of [T7] . 

3.2. Semiclassical expansions 

If A„ is deep inside the interval (—1, 1), the saddle points fl28|) shift from the real axis into 
the complex circle \z'\ = 1, and become two well separated complex conjugate numbers. 
Then n ^ 1 at /i^ — )■ -|-0, and one can easily obtain the semiclassical expansion for 
En{P) from the saddle-point asymptotics of the integral (!27|) . To the leading order in 

fj, ~ hz, we get 

E„(P) = J - 4A„ a cos(P/2), (32) 

A„ = - cos On, 

sm9n - 6'„cos6'„ = /ivr (n - 1/4) + arg(l ae*^") + 0{i/), 
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sm / dz [a + z ) exp 

\ J J-i 



0. (33) 



where parameters a and /i are given by (fT^ . 

On the other hand, if A„ is well above 1, integral in f l27|) being determined by the 
saddle point — A„ + a/A^ — 1, lying in the real interval (—1,0). However, there are 
still two saddle point contributions to the integral in (!27|) coming from the upper and 
lower edges of the contour C, see Fig. lAll Since these contribution differ only by the 
phase factors exp(^i7rA„//i), equation fl27|) can be written at /x — )■ +0 with exponential 
accuracy as 

-Anlog|2:| - { z - z'^)/2 

Equating to zero the first factor in (!33i) leads the equidistant Zeeman ladder 

A„ = /i(n + l), (34) 

EniP) = J + hin + 1), (35) 

with n = 1, 2, 3, . . . Corresponding bound states can be considered semiclassically as two 
well separated localized kinks moving back and forth without mutual collisions |16j . 

Putting to zero the integral in fl33|) gives the energy of the well localized kinetic 
mode (|TT1) modified by the longitudinal magnetic field /i^: 

1 + a'^ fi 

2a ^ a2 - 1 



± -r U, Lb Ll ^ / 9 

Xk^n = — + + O /i' 



EUP) = J +(/?'- 2/3 cos P) + '^rr^^osp + (36) 



(/3'-2/3cosP) 



2 



2^ 



(/?' - 2(3 cos P)2 - [2a cos(P/2)]2 + 
Fig. |3] shows the same spectra as in Fig. [21 together with six asymptotical curves. 
The filled area displays the region of the continuous spectrum at hz = 0. Five curves 
in the bottom display the spectra of five lightest modes calculated from the low-energy 
expansion fl29|) . This expansion does not hold outside the filled region. The curve 
crossing the Zeeman ladder in the top of Fig. [3] represents the kinetic bound-state mode 
determined from ( |36|) . It is clear, however, that this are in fact the avoided crossings 
with exponentially narrow gaps. 

Note, that all above asymptotical formulas cannot be used in the crossover region 
near the upper bound of the filled area. 



4. Dynamical correlation function 

Besides the magnetic excitation energy spectra, the inelastic neutron scattering allows 
one to measure the dynamical correlation function [T| [T7] 

oo 

s''{p,u) = syy{p,u) = J2\{'^n{P)\s'om's[u-E4P)], (37) 

n=l 

where |0) = | . . . ttt • • •) is the ferromagnetic ground state, Sj = (7^/2 is the x-spin 
operator at the site j, and the sum extends over all eigenstates of the Hamiltonian with 




Figure 3. The same spectra as in Fig. [2] witli added asymptotical curves explained 
in the text. 



the momentum P. Since the operator Sj inverts just one spin at the site j in the chain, 
it maps the ferromagnetic vacuum into the two-kink state, 

5J|0) ^ . . . ttttt . . .) = ^ I • • . ttitt |J, 1). 

Accordingly, the dynamical correlation function fl57|) for model (E]) takes the form 

S'%P, UJ) = -Y, In{P) S[u - E4P)], (38) 



4 

n=l 



n=l 



where In{P) is the relative intensity of the n-th mode, 

/ iP) - '^"(^ = ^'-P)'' (39) 

and ipni^, P) are the eigenfunctions of Hamiltonian ([7]). 

Equations ( 137|) -( 139|) are written in the assumption that Hamiltonian (|3]) has only 
the discrete spectrum. These relations can be easily modified in the case = 0, where 
the continuous spectrum exists: 

S^%P, CO) = \hin{P) 5[u - E,UP)] + ^ I{k, P) 6[u - E{k, P)], (40) 

4 271 Jq 
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where 

In ;f I „ I ^ 1 

(41) 
(42) 

1 + 2a cos + a2 ' ^ ^ 

and parameters a and 7 are given by f|T5|) and f lT2|) . respectively. Fig. H] displays the 





Figure 4. Intensity I{k,P) of the continuous modes at hz — in the (P, i?)-plane 
determined from (H^ . ([T^ . Darker regions correspond to larger intensity. 

density plot of the intensity I{k, P) of the continuous modes at hz = in the (P, E)- 
plane determined from (H2l) . (|T2l) for the values of the rest parameter given by ( l26l) . 
Formulae f l37|) - fl39|) can be directly used at nonzero longitudinal magnetic field 



hz > 0, since the spectrum is discrete in this case. Using the results of Appendix B 
the relative intensities of the discrete modes can be expressed in terms of the Bessel 
function 



(43) 



1- 


\ JuHlt^) 1 


Y 






_J^+l(l//i)_ 
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where /z is given by (IT^ . and z/„ is the n-th solution of equation ([22]) • 




P 



Figure 5. Intensities of the discrete modes In{P) at a weak longitudinal magnetic 
field according to (|43p , with the parameter values . The solid curve represents the 
intensity Ikin{P) of the kinetic mode at —0, according to (|1T|) . 

Figures [5] and [6] show intensities In{P) of discrete modes calculated from f H3|) . f|T5|) 
with the parameter values f l26|) . Intensities /n(0) decrease with increasing n. Intensities 
of 30 lightest modes are plotted versus the momentum P in Fig. El and those for 
n = 1,2,3,4,5,18 are shown in Fig. |6l The intensity of the kinetic mode Ikin{P) at 
/i^ = is shown in Fig. |5]by the solid curve. 

The obtained results are summarized in Fig. [71 It shows the same energy spectra 
as those in Fig. [2l The dispersion curves were calculated by use of equations (l24l)-(l26l). 
The darkness of the dispersion curves in Fig. [7| characterizes the relative intensities ( l43l) 
for the corresponding modes. This figure looks quite similar to Fig. 3B of ref. [1], which 
displays the experimentally observed neutron scattering spectra in CoNb206. 

5. Conclusions 

In this paper we obtain the exact solution of the phenomenological model, which was 
proposed by Coldea et al [T7] to describe the spectra of two-kink bound states in 
the entire Brillouin zone observed in the quasi-lD Ising ferromagnet CoNb206 in the 
inelastic neutron scattering experiments [1]. The model Hamiltonian acts in the space 
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Figure 6. Intensities In{P) for modes witli n = 1, 2, 3, 4, 5, 18 at nonzero h according 
to P5|) . ([T5|). Intensities /„(0) decrease with increase ti. 



of two-kink configurations of the spin chain, and parametrizes hoppings of kinks along 
the chain, their short range interaction, and the hnear long-range attracting potential 
leading to confinement of kink into pairs. We express the dispersion law En{P) of 
two-kink bound states and relative neutron scattering intensities In{P) by the magnetic 
modes in terms of the Bessel function. The preliminary analysis shows [H], that the 
obtained spectra En{P) plotted in Fig. [2] display a good quantitative agreement with 
the experimental data. It would be interesting to perform the detailed comparison of 
obtained exact solutions [both for energies -E„(P) and intensities /n(-P)] with the neutron 
scattering data in CoNb206. 

Let us note, that the dispersionless modes forming the almost equidistant ladder 
with energies ( 135|) and momenta near the zone boundary have very small intensities 
In{P), as it is seen from Fig. [71 The reason is quite simple. It is well known, that the 
linear potential causes localization of a sole kink in the discrete spin chain, similarly to 
localization of an electron in the isolated conducting zone by the uniform electric field 
[19]. The dispersionless modes ( 135|) correspond to large enough clusters of down-spins 
bounded by two well separated localized kinks, which can be hardly excited in the up- 
spin ground state by the neutrons scattering. Perhaps, such modes could be effectively 
excited by application of the oscillating magnetic field having the resonant frequency 
oj = hz/h. 
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Figure 7. Energy spectra of 30 lowest modes according to 
of the curves characterizes the intensities (|43p of the modes. 
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Appendix A. Solution of the eigenvalue problem 

In this Appendix we describe the exact solution of the eigenvalue problem f ll4p . which 
we rewrite as 



2 J ^ ' ' 2 
with the Dirichlet boundary conditions 

lim ^{l) = 0, 

i— >+oo 

ipiO) = 0. 



(A.l) 

(A.2) 
(A.3) 
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The following normalization condition for the eigenstate will be chosen: 

^(1) = -2. (A.4) 

Let us skip for a while the boundary condition ( lA.Sp and consider the generating 
function g{z) for ipil), 

oo 

1=1 



(A.5) 



This function should be analytical inside the circle \z\ < 1, subject to the boundary 
condition 



9(0) = 0, 

and satisfy the differential equation 

dg{z) 



Hz) - ^]9{z) + tiz- 



dz 



V{z) 



following directly from (1A.1I) . (IA.4I) . ( ]A.5I) . Here e{z) is given by ( IT9l) . and 



Viz) 



l + z 



a + 



The solution of equation (1A.7P satisfying (1A.6|) reads as 



where 



9[z) 



iHz) 



dz' 



fJ'Z' 



-Xhgz - -{z 



-Viz')exp\-[Hz')-Hz)] 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



(A.IO) 



and the branch of the logarithm is fixed by the condition J^{1) = 0. 

At positive /i, the integration path in the 2;'-plane in flA.Qp should approach the 
origin z' = from the left half-plane Kez' < to provide convergence of the integral. 
Figure lAll shows two allowed integration paths 71 and 72 in the right-hand side of (lA.OP 
for the case z > 0. Integration along each of them should give the same function g{z). 
This requirement leads to the following constraint (cf. equations (56), (57) in ref. 

dz' 



c 



— V{z')exp\-[Hz')]\=0 



where the integration path C = 72 — 71 is shown in Figure 
into (lA.lip gives the constant ip{0): 



^(0) = -2 



a + 



(A.ll) 
Substitution of (lA.SP 

(A.12) 



Ju+l{t^ ^) 

where p = —A///, and Jiy{x) is the Bessel function. We have taken into account the 
integral representation of the Bessel function 



Joii^x) 



f ^z"eMx{z''-z)/2], 

Jc ^TTl Z 
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which reduces to the well known form (see formula (5) in page 15 in [20j) after the 
change of the integration variable z = 1/u. 

Application of the Dirichlet boundary condition (IA.3I1 leads then to the equation 



J,(/i-^) + a Ju+i{fi-^) = 0, (A.13) 

which solutions z/„ determine the spectrum A„ = —jJifn of the problem ( 1A.1[) -( |AI3!) . 

Fig. IA2I shows the plot of the left-hand side of equation ( ]A.13I) versus v ai a = —0.5, 
/X = 0.05. 

Appendix B. Relative intensities of discrete modes at /i > 

In this Appendix we prove the following 
Statement. 

Let the set of real numbers •^)}/So solves the linear problem 

, + ^, + ^^^„,,)_*(i±M±jW^ = 0. i = l,2...., (B.l) 



2 

with the boundary condition 



lim i){l,X) = 0. (B.2) 
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Figure A2. Plot of the function J,y(/i ^)+a Ju+i{n ^) versus v &t a = —0.5, ^ = 0.05. 



Coefficients a, X, are supposed to be real, and /i > 0. 

Let us define the intensity /(A) corresponding to this solution as 

[^(1,A)]^ 



/(A) 



Then 



Proof. 



J(A) = 2/i 



-1 



(B.3) 



(B.4) 



u=—X/fi 



Since definition ( IB. 31) of /(A) does not depend on the normalization of the set 
A)}^Q, we shall fix the latter by the condition 

7/>(l,A) = -2, (B.5) 

without loss of generality. Consider the solution {V'(^,A')}^q of the problem (IB.ip . 
(IB. 21) . (IB. 5 1) in which A is replaced by A'. In particular, instead of (IB. II) we get 



-X' + IJ,l + 



/ = 1,2,... 



(B.6) 



Let us multiply equation ( ]B.1|) by ipi^l, A'), and equation ( ]B.6P by ■?/'(/, A), then subtract 
one equation from the another and sum the result over all natural /. As the result, we 
obtain 

(A' - A) A) ^(/, A') = - [V(l, A')^(0, A) - ^(1, A)V^(0, A')] , (B.7) 



1=1 



or 



^V^(/,A)^(/,AO 



1=1 



1 

2(A' - A) 



:i,A'MO,A)-^(l,AMO,A0]. (B.8) 
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Taking into account the chosen normahzation condition iIj{1,X) = iIj{1,X') 
proceeding in (IB.8|) to the hmit A' — )■ A, we get 



1=1 



rf^(0,A) 
dX ' 



Combining f lR9|) with f[02|) yields 



1=1 



-2, and 



(B.9) 



(B.IO) 



Substitution of f lR5|) and flRTOj) into flR3|) leads finally to the result flR4|l . 
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